In [1]:
# Standard setup block for running Python code
from __future__ import division, print_function
import os
if os.path.split(os.getcwd())[-1] == "Lab notebooks":
os.chdir("../../")
print("Moved to experiment root directory")
from Modules.processing import *
from Modules.plotting import *
plt.style.use("Config/plotstyle.mplstyle")
%matplotlib inline
Got to the lab about 2:30 PM.
Homed turbine, tow, y-, and z-axes.
Tank depth is 2.44 m.
Mixed up seeding from surface with squeegee like yesterday.
2:48 PM -- Did two dummy tows to mix up tank further.
2:53 PM -- Restarting Wake-1.0-0.625 at run 16.
5:32 PM -- Finished Wake-1.0-0.625. Backed up data. z-axis at lower limit RPOS is -0.03600 m. y-axis is -1.52566.
5:37 PM -- Re-homed y- and z-axes. Starting Wake-1.0-0.75
9:09 PM -- Stopping Wake-1.0-0.75 at run 35. At negative limits RPOS(z) = -0.0359869
and RPOS(y) = -1.52510
. Backing up data.
In [77]:
Section("Wake-1.0-0.75").process(nproc=1, nruns="new")
In [78]:
plot_meancontquiv()
In [71]:
w = WakeProfile(1.0, 0.75)
w.plot("mean_u", "-o", preliminary=True)
In [5]:
r = Run("Wake-1.0-0.625", 29)
plt.plot(r.time_vec, r.u)
plt.show()
In [6]:
length = 0.5
diameter = 2.9*0.0254 # inches to m
area = length*diameter
cd = 1.2
r_force = 0.25
rho = 1000.0
u = 1.0
tsr = 5.5
f = 0.5*rho*area*cd*(tsr/2*u)**2
torque = f*r_force
print(f, torque)
In [7]:
Run("Perf-1.2", 5).torque_all.max()
Out[7]:
In [70]:
# Setting things up
setup_time_sec = 3600
n_setup_times = 4
# Seconds per tow at each speed
seconds = {0.8 : 90 + 200,
1.0 : 84 + 240,
1.2 : 80 + 300}
# Wake runs
runs_per_wake = 43
nwake = (43 - 33)/43
wake_speed = 1.0
# Runs for strut torque and strut covers
strut_covers_runs = 22
strut_torque_runs = strut_covers_runs*2
# Tare drag and torque
tare_drag_runs = 14
tare_torque_runs = 30
tare_drag_dur = 200
tare_torque_dur = 120
total_time_sec = n_setup_times*setup_time_sec + runs_per_wake*nwake*seconds[wake_speed] \
+ strut_covers_runs*seconds[wake_speed] + strut_torque_runs*tare_torque_dur \
+ tare_drag_runs*tare_drag_dur + tare_torque_runs*tare_torque_dur
total_hours = total_time_sec/3600
hours_per_day = 10
total_days = total_hours/hours_per_day
print(total_hours, total_days)